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Abstract. Absolute proper motions, their corresponding errors and membership probabilities of 501 stars in the 
intermediate-age open cluster NGC 2548 region are determined from MAMA** measurements of 10 photographic 
plates. The plates have the maximum epoch difference of 82 years and they were taken with the double astrograph 
at Z6-Se station of Shanghai Observatory, which has an aperture of 40 cm and a plate scale of 30" mm"'^. The 
average proper motion precision is 1.18 mas yr~^. These proper motions are used to determine the membership 
probabilities of stars in the region. The number of stars with membership probabilities higher than 0.7 is 165. 
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1. Introduction 

The open cluster NGC 2548, also known as M 48, in 
Hydra (a2ooo=8''13'°48", ^2000 = -5°48') seems to be an 
intermediate-age open cluster with an estimated distance 



of 630 pc (Pesch |l96l| ) or 530 pc (Claria |1985|). It has 
I not been the object of any recent complete astrometric 
or photometric study (Ebbinghausen 1939| ; Li 1954), in 
spite of being an extended object with an apparent dia- 
meter of 30' (Trumpler |1930| ) or even 54' (ColHnder [1931D 
and brilliant enough to b e in t he Messier list (s. XVIII) 
with number 48 (Messier 1850 ). It was even considered 
inexistent for several years owing to the fact that Messier 
mistook its real coordinates owing to a change of sign in 
its relative position with a final difference in declination 
of 5°. 



There is no feasible estimation of its age but it seems 
to be an intermediate-age open cluster, around logt = 8.5 
(Lynga 1987 ), with a slightly poorer CN abundance than 
the giants of the Hyades but significantly r icher than the 
K gia nts o f the solar neig hborhood (Claria |1985| ; Twarog 
et al. |19971). 
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* Table 5 is only available in electronic form from CDS 
via anonymous ftp 130.79.128.5, and by e-mail request to 
Lola.Balaguer@am.ub.es 

** MAMA ( Machine Automatique a Mesurer pour 
I'Astronomie, tittp://dsmama.obspm.fr) is developed and 



operated by INSU/CNRS and Observatoire de Paris. 



Proper motions of bright stars in this region were pub- 
lished by Ebbinghausen ( 1939 ) with four pairs of pl ates o f 
only a maximum epoch difference of 28 years. Li (1954) 
published positions and proper motions in the field of 
NGC 2548 based on three plates taken with the 40 cm 
astrograph at Shanghai Zo-Se station with a maximum 
epoch difference of only 14 years. 

In this paper we determine, for the first time, precise 
absolute proper motions of 501 stars within a 1?6 x 1?6 
area in the NGC 2548 region, from automatic MAMA 
measurements of 10 plates, five of which are newly taken. 
The estimated membership probabilities lead us to a com- 
plete astrometric study of the cluster area. Section 2 des- 
cribes the plate material and its measurement as well as 
the proper motion reduction procedure and results with 
comparisons with the Hipparcos and Tycho-2 catalogues. 
Section 3 accounts for the membership determination. 
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Section 4 is devoted to the analysis of results. Finally, 
a summary is presented in Section 5. 

2. Plate measurement and proper motion 
reduction 

2.1. Plate material and measurements 

Ten plates of the NGC 2548 region were available. They 
were taken with the double astrograph at the Z6-Se station 
of Shanghai Observatory. This telescope, built by Gaultier 
in Paris at the beginning of the last century, has an aper- 
ture of 40 cm, a focal length of 6.9 m and hence a plate 
scale of 30" mm^^. The size of the old plates is 24 cm by 
30 cm, or 2?0 x 2?5, and that of the new ones is 20 cm by 
20 cm, or 1?65 x 1?65. 

The oldest plate was taken in 1916, and the newest 
ones in 1998. Relevant information on these plates is 
shown in Table 0. The hour angles are not provided in 
Table |^ because the starting exposure time of the old 
plates was not recorded. 

All plates were measured at the Centre d' Analyse des 
Images at the Observatoire de Paris, using the high preci- 
sion microdensitometer " Machine Automatique a Mesurer 
pour r Astronomic" (MAMA). This device has a super b 
optical and mechanical performance (Guibert et al. |1990|) . 
It uses a quartz-iodine illuminating source, whose trans- 
mitted light is detected by a reticon, 1024 pixels large, 
with a pixel size of 10 /im, and the absolute accuracy of 
the measurements is 0.6 /im (Soubiran 1992 ). After the 
plates were scanned, the resulting images were stored in 
a grid of 19 x 19 sub-images for each plate. Once every 
plate was digitized, we identified all point sources in these 
361 frames. The source extraction was performed on each 
frame using SExtractor (Berlin & Arnouts 1996), a soft- 
ware dedicated to the automatic analysis of astronomical 
images using a multi-threshold algorithm allowing good 
object deblending. To improve the accuracy of the mea- 
surements, we chose to retain only rather bright objects 
with a signal-to-noise ratio of at least 12. 

During the scanning, MAMA includes in the catalogue 
not only real images, but also the plate grid, emulsion 
flaws, plate annotations and scratches. A first step is re- 
quired to reject spurious detections, most of which can 
be recognized with analysis of the object shape, by re- 
moving the plate grid, annotations and big scratches by 
visual comparison with the original plates. The shape pa- 
rameters given by SExtractor can be used to clarify the 
remaining stars but the most direct way to reject spurious 
detection is by comparison between plates, searching each 
object from one plate to the other, and retaining only the 
paired data. There is a total of 596 stars measured with 
a limiting magnitude Bt around 14. This limiting mag- 
nitude was roughly estimated from the stars in common 
with Tycho-2 catalogue. 

The detection of 182 PPM and ACT stars on the whole 
plate allowed us to determine the star positions in a com- 
mon system. Thus for each scanned plate, we obtained a 
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Fig. 1. The number of stars vs the number of available 
plates 



preliminary astrometric catalogue to perform initial cross- 
identification between plates. 



2.2. Proper motions 

The absolute proper motions for 596 stars in the region 
of NGC 2548 were reduced on the basis of the MAMA 
measurements following the central overlapping procedure 



(Russel |1976| Wang et al. |1995| , |1996| , |2000[ ). This method 
determines the plate-to-plate transformation parameters, 
the star motions and their errors simultaneously. Stellar 
positions and absolute proper motions were reduced from 
a catalogue used as the original data for the first itera- 
tion. As initial catalogue, 265 stars from the Tycho-2 ca- 
talogue at epoch J2000 (H0g et al. 2000) were chosen on 
the basis of the results of the PPM and ACT astrome- 
tric catalogue given by MAMA. To select the best plate 
constant model, we used Eichhorn & Williams' criterion 



(Eichhorn & Williams |1963[ , Wang et al. |982|) and ob- 
tained a model with six linear constants on coordinates, 
a magnitude and a coma term, and a magnitude distor- 
tion term. Magnitudes used were the instrumental magni- 
tudes. All the proper motions are constrained by having at 
least one measurement from the modern epoch plates, i.e. 
taken in 1998. The whole process is iterated until the resul- 
ting proper motions converge. We iterate the process until 
mean differences in position are smaller than 1.1 mas, the 
r.m.s. smaller than 3.6 mas and the differences in proper 
motion below 0.1 mas yr^^, yielding a final outcome of 
501 stars. 

Table |^ shows mean precisions of final proper mo- 
tions for stars in the NGC 2548 region detected on dif- 
ferent numbers of measured plates (greater than 2). The 
units of the proper motions and their precisions all along 
arc mas yr^^. The precision of the final proper motions 
strongly depends on the number of plates. Figure ^ gives 
the number of stars for which various numbers of plates 
are available. More than 90% of proper motions were ob- 
tained from at least 5 plates. 
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Table 1. Plate material 



Plato 


Epoch 


Exp. Time 


Plate center 


Plate 


size 


11. of 


Id. 


(1900+) 


mill 


(aj2000 


5j20oo) 


cm 


stars 


CL422 


16.03.24 


90 


8!'2283 


— 5? 779 


24 X 


30 


450 


CL534 


30.03.21 


90 


oh Q/; 


CO '7'7C 
— O. MO 


24 X 


30 


556 


CL535 


30.03.28 


90 


8?'2291 


-5?767 


24x 


30 


566 


CL56006 


56.03.14 


60 


8!'2288 


-5? 779 


24 X 


30 


577 


CL56007 


56.03.16 


60 


8!'2372 


-5? 779 


24x 


30 


558 


CL98004 


98.04.03 


30 


8?'2447 


-5?958 


24 X 


30 


548 


CL98047 


98.12.16 


30 


8?'2227 


-5?796 


20x 


20 


268 


CL98Tian 


98.12.25 


30 


8?'2255 


-5?804 


2Qx 


20 


268 


CL98Chon 


98.12.25 


30 


8?'2255 


-5?804 


20x 


20 


298 


CL98GU 


98.12.16 


30 


8!'2227 


-5?796 


2Qx 


20 


432 



Table 2. Mean precisions of proper motions as a function of the number of plates in the NGC 2548 region. (Units are 
mas yr~^.) N is the number of stars 



plates 


iV 


^fla COS S 






3 


4 


5.29 


2.90 


6.30 


4 


8 


2.92 


3.12 


4.44 


5 


46 


1.76 


1.11 


2.19 


6 


119 


1.19 


0.89 


1.53 


7 


27 


1.15 


0.76 


1.42 


8 


20 


0.90 


0.70 


1.17 


9 


50 


0.74 


0.55 


0.94 


10 


199 


0.59 


0.46 


0.77 



The r.m.s. errors on proper motions for more than 90% 



of stars are e. 



, cos S 



5= 0.92 mas yr 



and 1.18 mas yr" 



where — 



0.68 mas yr ^ 
In 



i.iu liicxo , „ii^i^ — Y V„ COS 5 ^ "^Mci 

the most precise case, the r.m.s. errors are 0.77 mas yr^^ 
for stars with 10 plates (40% of stars). 

Figure |^ shows the distribution of r.m.s. proper mo- 
tion errors with the number of stars: N versus e^^cosS, 
e^^ and e^. Thus the precisions of the proper motions of 
stars in the region of NGC 2548 obtained in this study are 
relatively high, thanks to the quality of the stellar images 
taken with the 40 cm double astrograph and the excellent 
positioning behavior of the MAMA scanning machine. 

Figure ^ gives fj,a cos 6, fxs and their errors with respect 
to instrumental magnitude. 

Our absolute proper motions and their errors are com- 
pared with those of Tycho-2 catalogue in Fig. ^. We 
would like to highlight the precision of the proper mo- 
tions derived in this paper. Mean differences in the sense 
ours minus Tycho-2 are -0.123 (cr = 2.112) and -0.203 
(a = 2.158) in fia cosS and fj,s, respectively. No apparent 
systematic residuals were found as a function of magni- 
tude. A linear fit to the proper motion data gives us: 

fia cosS = 0.187(±0.133)-l- 0.982(±0.009) ■ {^a cos^)Tyc2: 
r = 0.988 

fis = 0.201 (±0.132) + 1.016 (±0.008) • {fis)Tyc2 ; 
r = 0.992 

being r the correlation coeficient. 

Only 8 stars were found in this region from the 





10 11 12 13 14 




10 11 12 13 14 



Hipparcos catalogue (ESA 1997). By comparing our abso- 
lute proper motions with these common stars, the mean 
differences are (in the sense ours minus Hipparcos) —0.610 



Fig. 3. Proper motions (top) and their errors (bottom) vs 
instrumental magnitude. Null errors are from proper mo- 
tions calculated with only two plates. (Units are mas yr^^) 



(cr = 1.710) and —0.198 (tr — 1.658) in fiaCOsS and fis, 
respectively. We obtain the following linear fit: 

fia cos(5 = -0.458(±0.756)+ 1.005(±0.013)- (/i^ cos(5)hip; 
r = 0.999 

= -0.470 (±0.694) + 1.011 (±0.014) • {fis)mp ; 
r = 0.999 
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Fig. 2. The number of stars vs the r.m.s. errors on proper motions in mas yr ^ 
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Fig. 4. Proper motions and their errors from this pa- 
per compared to those in Tycho-2 catalogue. (Units are 
mas yr"'^) 



3. Membership determination 

Accurate membership determination is essential for fur- 
ther astrophysical studies of clusters. The fundamental 



mathematical model set up by Vasilevskis et al. (1958) 
and the technique based upon the maximum likelihood 



We used a maximum likelihood method with a 9- 
parametric Gaussian model for the frequency function, as 
follows: 

$ = $c + = ■ 0c + n/ ■ 

where (pc and (pf, are the probability density functions of 
cluster members and field stars respectively in the proper 
motion space, with Uc the normalized number of cluster 
stars, and nj the normalized number of field stars. 

The probability density function for the i-th star of 
the cluster can be written as follows: 



S) 



1 



2^(a2 + e2 



exp ■ 



, )l/2(a2 +£2 , )i/ 

[{fla COS d)i - {jJLa COS 5)cY 



-(Mo 



3 5). 



principle developed by Sanders (1971) have since been 
continuously refined. 

An improved method for membership determination 
of stellar clusters based on proper motions with different 



where [(/iaCosJ)^ ,(/i5)i] are the proper motions of the 
i-th star, [(/ia cos (5)c, (/i5)c] the cluster mean proper 
motion, CTc the intrinsic proper motion dispersions of 
member stars and [e(^^ cosA')^ ^^(/^s)^] the observed errors 
on the proper motion components of the i-th star. 

Analogously, for the field we have: 

^ ' ■ (Mc<COBi5)j: (fI„COBl5)i' ^ ( ,5 ) / ( f» « ) i 



observed precisions was developed by Stetson (198C) and 
Zhao & He ( |1990D . Zhao & Shao ( |1994| ) then added the 
correlation coefficient of the field star distribution to the 
set of parameters describing their distribution on the sky. 



exp 1 


1 




COS (5)^ —(l-la COS 5) f]'^ 


2(1-72) 




a cos 5) j.- (l-ia cos 5)j 


27[(Aia cos 5)i- 


- (/icv COS 5) j 






i (flQ COS l5) J (nc 







This model has been frequently used (Wang et al. 1995 
199^, |2000|). 



where [{pa cos 6) f , {fj,s) f] are the field mean proper mo- 
tion, [<7(^^ cos(5)j> 7 f (/^s)^] the field intrinsic proper motion 
dispersions and 7 is the correlation coefficient. 
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4. Results and discussion 

The unknown parameters for the assumed $ distribu- 
tion are [ric, (/^q, cos(5)c, (/^(5)c, fc] for the cluster and 



[(/XaC0S(5)/, {fJ.s)f, cr. 



(Ha cos (5)/ J '^{fJ.s)fj 



7] for the field popu- 



lation. Membership probability of the i-th star belonging 
to the cluster can be calculated from 



By applying the standard maximum likelihood 
method, we obtained the 9 distribution parameters and 
their corresponding uncertainties, shown in Table ^, where 
the units of the proper motions and proper motion dis- 
persions are mas yr~^. The quality of the fit should be 
optimized near the proper motion region occupied by the 
cluster stars, where the model is most crucial to providing 
reliable membership determination. Outliers cause the es- 
timated distribution of field stars to be flatter than the 
actual one, thus increasing the final probability of mem- 
bership (Kozhurina-Platais et al. 1995, Cabrera-Cafio & 
Alfaro |1984 Zhao & Tian |1985| , Zhao et al. ^982|). To 
minimize the effect of high proper-motion field stars in 
the model, we restricted the membership determination 
to the range \fj,\ < 30 mas yr~^. 

Table 5 lists the results for all 501 stars in the region 
of the open cluster: column 1 is the ordinal star number; 
columns 2 and 3 are aj2ooo and (5j2ooo; columns 4 and 6 
are the respective absolute proper motions (/Xq. cosJ, /i^); 
columns 5 and 7 are the standard errors of the proper 
motions; column 8 is the membership probability of stars 
belonging to NGC 2548 (Pc); column 9 is the instrumen- 
tal magnitude given by SExtractor; column 10 gives the 
number of plates used and column 11 the identification 
number in Tycho-2 for the stars in common. 

The cross-identification of our 8 common stars with the 
Hipparcos catalogue is shown in Table ^. Comparison with 
the membership determination calculated by Baumgardt 



et al. ( 2000) for the seven stars in common shows complete 
agreement. 

Figures || and || show the proper motion vector-point 
diagram and the position distribution on the sky for all the 
measured stars respectively, where " o " denotes a mem- 
ber of NGC 2548 with Pc > 0.7, and all other stars are 
considered field stars indicated by . 

The cluster membership probability histogram (Fig. ^ 
shows a very clear separation between cluster members 
and field stars. The number of stars with membership 
probabilities higher than 0.7 for NGC 2548 is 165. The 
average cluster membership probability is 0.90, giving a 
contamination by field stars not larger than 10%. If these 
stars were considered as belonging to the field, the ap- 
parent deficiency of field stars in the central region would 
vanish. Likewise, the contamination by cluster stars in the 
field is not larger than 6%. 

The present results for this cluster should be comple- 
mented with photometry and radial velocity studies. To 
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Fig. 7. The histogram of cluster membership probability 
of NGC 2548 



this end, our group is working on a complete photometric 
study of the region. 

4.1. Effectiveness of membership determination 

Contamination by background and foreground objects due 
to the observational projection effect cannot be avoided. 
Following Shao & Zhao (1996), we can assess the effec- 



tiveness of our membership determination, which is set 



^ELm(^) [l-Pc(^)]} 

Eti^c(*)Ef=i [i-Pc{i)] 



E = 1 



where the greater the more effective the membership 
determination. 

We thus obtained an effectiveness of membership de- 
termination of 0.77 for NGC 2548. Figure 3 of Shao & 
Zhao ( 1996| ) paper shows that the effectiveness of member- 
ship determination of 43 open clusters ranges from 0.20 to 
0.90 and the peak valu e is .55. Compar ed wit h previous 
studies (Sha o fc Zhao 199£; Tian et al. 1998 ; Balaguer- 



Niifiez et al. 1998), the effectiveness of membership deter- 
mination for this open cluster is significantly high. 

4.2. Space velocity and galactic orbit 

By combining our absolute proper motion results with ra- 
dial velocity (Piatti et al. 1995D , and the age and distance 
of the cluster from the Catalogue of Lynga (1987), we 
determined the space velocity and the galactic orbit of 
NGC 2548. 

To obtain the velocity of the cluster in the galacto- 
centric frame, we assumed the motion of the Sun in the 
LSR to be ([/,!/, 1^)0 = (9.7,5.2,6.7) km s"! (Bienayme 
1999). We adopted the current lAU standard values of 
Vlsr = 220 km s~^ for the local circular rotation velocity 
and i?0 — 8.5 kpc for the Sun galactocentric distance. The 
resulting galactocentric position and velocity of the clus- 
ter are {x, y, z) = (-8.89, -0.44, +0.16) kpc and [U, V, W) 
= (0,222,7) km s^^. These vectors, together with the 
galactic gravitational potential model of Allen & Santillan 
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Table 3. Distribution parameters and their uncertainties for NGC 2548. The units of /i and a are mas yr~^ 







/Iq cos (5 






o" n„ cos & cr^g 7 


NGC 2548 


0.382 


-1.41 


1.64 


1.23 






±0.025 


±0.12 


±0.12 


±0.08 




field 




-4.89 


-1.63 




7.37 8.21 -0.28 






±0.07 


±0.57 




±0.04 ±0.29 0.03 



Table 4. The cross-identification of stars in common with the Hipparcos catalogue and the membership determination 



by Baumgardt et al. (2000). Comparison with the results from this paper (Pc) gives complete agreement 



Table 5 


Hipparcos 


Tycho-2 


BDA 


Member 


Pc 


257 


40110 


4859.00078.1 


366 


N 


0.00 


140 


40238 


4859.00250.1 


1005 


M 


0.87 


139 


40254 


4859.00036.1 


1073 


M 


0.77 


336 


40281 


4859.00921.1 






0.00 


133 


40302 


4855.01706.1 


1320 


M 


0.78 


234 


40348 


4859.01156.1 


1560 


M 


0.88 


162 


40362 


4859.00092.1 


1628 


M 


0.84 


42 


40498 


4856.00072.1 


2184 


N 


0.00 



40 



30 



"3 



20 
10 


-10 
-20 
-30 



-40 



-40 -30 



-20 -10 10 

|i„cos5 



20 30 



40 



Fig. 5. The proper motion vector-point diagram of NGC 2548. (Units are mas yr ^.) ("o" denotes a member of 
NGC 2548 with > 0.7, "±" a field star with Pc < 0.7) 
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Fig. 6. The position distribution of stars in NGC 2548 area, ("o" denotes a member of NGC 2548 with > 0.7, "+" 
a field star with Pc < 0.7) 



( |l99lD with the estimated cluster age of logt = 8.5, de- 
termine the orbit of the cluster in the Galaxy. The re- 
sulting orbit is characterized by an eccentricity of 0.022, 
radial oscillations between distances from the galactic cen- 
ter of 8.67 and 9.04 kpc and vertical oscillations between 
—0.19 and 0.19 kpc. The cluster has made 1.18 revolutions 
around the galactic center and 8 crossings of the galactic 
disk in its lifetime. Assuming the distance of Claria (1985) 
does not affect the results. We currently observe the clus- 
ter near its maximum distance from the plane. This is 
consistent with the expectation of a higher probability of 
finding the cluster near its maximum \z\, where the cluster 
will spend more time due to its vertical motion. 

5. Summary 

Proper motions and their corresponding errors of 501 stars 
within a 1?6 x 1?6 area in the NGC 2548 region were 
determined from automatic MAMA measurements of 10 
plates, five of which were newly taken. Comparison with 
Tycho-2 and Hipparcos Catalogues shows good agreement 
and underlines the precision of the proper motions de- 



rived in this paper. These proper motions are then used 
to determine membership probabilities of the stars in the 
region. We obtained 165 stars with membership probabi- 
lities higher than 0.7. 
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